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We review the recent progress in new lattice fermion formulations. We focus on the follow- 
ing three types which have possibility of improving lattice simulations. (1) Flavored-mass 
fermions are a generalization of Wilson fermions with species-splitting mass terms. In particular, 
staggered- Wilson fermions initiated by Adams have possibilities of reducing numerical costs in 
overlap fermions and the influence of taste-breaking in staggered fermions. (2) Central-branch 
Wilson fermions, in which additive mass renormalization is forbidden by extra axial symme- 
try, could enable us to perform Wilson-fermion lattice QCD without fine-tuning. (3) Minimally 
doubled fermions, which reduce the number of species by species-dependent chemical potential 
terms, realizes a ultra-local chiral fermion at the price of hypercubic symmetry. These setups 
reveal unknown aspects of lattice fermions, and we obtain a deeper understanding of lattice field 
theory. 
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1. Introduction 

Non-perturbative lattice calculations have contributed to a broad range of topics in particle and 
nuclear physics. It becomes more and more important in the era of new physics in LHC. However, 
practical problems in lattice fermions still prevent us from performing simulations efficiently. By 
now several fermion constructions have been developed, although all of them have their individual 
shortcomings. To perform more efficient and precise lattice simulations and to make progress in 
the frontier of high energy physics, it is of primary importance to develop lattice fermions with less 
numerical costs, less discretization errors, less tastes and better chiral property. Although effective 
ways of improving have been proposed and are used widely in the calculations, an interesting goal 
is to construct new fermion formulations free from the shortcomings. The research in this course 
will also give us a new insight into lattice field theory and feedback for the existing simulations. 

In this paper we review the recent progress on new fermion setups and discuss their appli- 
cations in lattice simulations. In Sec. ^ we discuss generalization of Wilson fermions by using 
species-dependent mass terms, which we call "flavored-mass" ||, ||]. In this study we will find 
new versions of Wilson, domain-wall and overlap fermions, some of which successfully have bet- 
ter dispersion relation and better eigenvalue spectrum [Q]. In particular staggered-Wilson (overlap) 
fermions Jj], |5]] have the possibility of reducing overlap computational costs [^, ^ and reducing 
staggered taste breaking effectively [§, ^]. We secondly discuss a different use of Wilson-type 
fermions. In Sec. || we will show that the Wilson fermion without non-hopping terms has an extra 
U(l ) symmetry, which forbids the additive mass renormalization [J7| [R]]. This axial symmetry is 
spontaneously broken with the associated NG boson emerging. By combining this idea to the fla- 
vored mass, we in principle obtain a two-flavor central-branch fermion which could suit two-flavor 
lattice QCD. In Sec. g we discuss another way of keeping chiral symmetry and reducing the number 
of flavors in the lattice fermion. "Minimal-doubling setup" reduces the species to two with keeping 



one exact chiral symmetry at the price of hypercubic symmetry [11, 12, 13]. To perform lattice 
QCD with this fermion, we need to tune the several parameters to restore Lorentz invariance in the 
continuum limit [ 14, |TJ]. Recently JI^ ] introduced a new view that this setup can be inteipreted as 



a lattice fermion with species-dependent chemical potential terms. Application of this fermion to 
in-medium QCD is discussed in [17]. In Sec. S we discuss the perspective on these formulations. 



2. Flavored mass and generalized Wilson fermions 

In this section we discuss generalizations of the Wilson fermion from the viewpoint of species- 
splitting mass terms, which we call "flavored mass terms". It leads to new types of domain-wall 
and overlap fermions. This course is initiated by Adams in [jl|], and the intensive researches have 
been done by now. This topic is roughly classified into two parts: generalization of Wilson fermion 
based on naive fermions and its generalization based on staggered fermions. We begin with the 
former one and move to the latter. 

2.1 Flavored mass for naive fermions 

In this subsection we introduce flavored-mass terms for naive fermions as generalizations of 
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Figure 1: Free Wilson Dirac spectrum. The degenerate spectrum of 16 species in naive fermions are split 
into five branches with 1, 4, 6, 4 and 1 flavors. 



the Wilson term |g|]. The Wilson fermion action is given by, 

n,fl n n,n 



(2.1) 



where = {T +}1 - r_ M )/2, C M = {T +pL + r_ M )/2 with T ± ^ n = £/„,± M y„± M . The free Dirac s P ec " 
trum for the Wilson fermion is schematically depicted in Fig. [I]. The degeneracy of 16 modes in 
naive fermions is lifted into 5 branches, to which 1, 4, 6, 4 and 1 flavors correspond. We emphasize 
the three important properties of the Wilson fermion, 75-hermiticity, hypercubic symmetry and the 
lattice laplacian form ~ a J dx 4 \j/ x AY x + 0(a 2 ). These can be criterions for generalization. 

Now we briefly look into flavor-chiral symmetry of naive and Wilson fermions by following 



[10]. As well-known the massless naive action possesses U(4) x U (4), which is regarded as rem- 
nant of the continuum flavor-chiral symmetry group for 16 flavors. This U{4) x U (4) is given 
by 



Yn — > ex P 



Here, ri +) and F 
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are site-dependent 4x4 matrices: 
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(2.2) 

(2.3) 
(2.4) 



where n 



1" 



n^+n v and v = Ep^u v n p ■ See [JlOp for details. Quark condensate 



or quark mass break this U(4) xU (4) down to the U (4) vector subgroup . We call T^~' as 
vector-type group and T^~' x as axial-type group. In the presence of the Wilson term this U (4) x 



U(4) invariance is broken down to the U(l) invariance under I4 in Eq.(2.3). This generator is 
vector-type, which means that the Wilson fermion loses all the axial(chiral) symmetry. 

Now we go on to the main theme "flavored-mass terms". In [|3j], it was shown that there are 
four nontrivial types of flavored masses for naive fermions, which satisfy 75-hermiticity, possess 
the hypercubic symmetry and becomes covariant laplacian with proper mass shifts. The four types 
are classified based on the number of transporters, where we name the 1-link case as vector (V), 
2-link as tensor (T), 3-link as axial-vector (A) and 4-link as pseudo-scalar (P), 

4 

My = £C M , M T = £ £ C M C V , M A = £ £ ]JC V , M P = £ n C f > < 2 - 5 ) 

[1 perm. sym. perm. sym. V sym.\i=\ 
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Figure 2: (Left) Dirac spectrum for Mp [|j], where we have two branches with 8 and 8 flavors. (Center) 
Dirac spectrum for Mj, where we have three branches with 6, 8 and 2 flavors. (Right) Dirac spectrum for 
My +Mj +Ma +Mp [ph, where we have two branches with 1 and 15 flavors. 



where Y*perm. means summation over permutations of the space-time indices. Y*perm. an( ^ Hsym. we 
defined as containing factors, for example 1/4! for Mp. The vector type My is the usual Wilson term 
up to the mass shift, and the Wilson term in ( |2.1[ ) is just given as Wn(4 — My)\j/ n . In the momen- 
tum space, they are transformed into M v — > cos , Mj — > cos cos p v , Ma — > cos p^ cos p v cos p p 
and Mp — > cos pi cos /?2COS p^ cos ^4. As shown in [g], these flavored masses are also expressed in 
the spin-flavor representation as Mp ~ (1 <g> (T3 <8> T3 (8> T3 (2) T3)). 

By introducing the flavored-mass terms % x M F \\f n (F=V, T, A, P) or their combinations into 
the naive action instead of the usual My, we obtain generalized versions of Wilson fermions. The 
species-splitting depends on the types or combinations. Mp splits the 16 degenerate modes into two 
branches with 8 and 8 flavors. The Dirac spectrum for DNaive +Mp is depicted in the left figure of 
Fig. ||. Mr splits them into three branches with 6, 8 and 2 flavors as depicted in the center figure 
of Fig. |[ Ma has the same species-splitting as the case of My (Wilson). Combinations of these 
flavored-mass terms are also useful to split the species in a desirable form. For example, the total 
sum of the flavored mass terms My +Mt + Ma + Mp splits the 16 species into two branches with 
1 and 15 flavors as shown in the right figure of Fig. ||. We note the hypercubic symmetry remains 



as long as we consider a combination of the forms in ( |2.5[ ). The remaining flavor symmetry also 
depends on the types. For example, Mp with the mass shift breaks the symmetries of naive fermions 
in Eqs.( |2.3[ )(2.4) into the subgroup 



r (+) 



14, (-1) 



»l+...+«4 



75, (-1) 



(2.6) 



Although all the symmetries in F^' are broken as with the usual Wilson, but the remaining group 
is larger. This symmetry group is interpreted as the subgroup of eight-flavor symmetry in the 
branches. What we showed here are just examples. The notion of flavored-mass gives us a wide 
class of Wilson cousins with desirable numbers of flavors. It is obvious that from the flavored-mass 
we can also obtain generalized domain-wall and overlap fermions. 

Brillouin fermion : In Ref.[Q], a way of improving dispersion relation in Wilson fermion was 
discussed independently. This is called "Brillouin fermion". Here the Wilson term is replaced by 
the other form of the covariant laplacian, which minimizes the breaking of the rotational symmetry 
near the Brillouin boundary. From the viewpoint of the flavored-mass terms, the Brillouin fermion 
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corresponds to the case of the total sum of the four flavored-mass terms as 

My + Mj + Ma + Mp, (2.7) 
which split the 16 species into 1- and 15-flavor branches. (We note that a similar form was referred 



to in the context of improvement of overlap [18]) We can utilize the 1 -flavor branch as Wilson 



fermion or overlap kernel. |Q] argues that the merits of this fermion is that not only the dispersion 
relation becomes more continuum-like, but the Dirac spectrum becomes more Ginsparg-Wilson- 
like. The calculation in the paper also shows that this setup is more effective as the quark mass 
is heavier. It is an interesting course to consider its application to the heavy quark systems. One 
concern about this formulation is that the four flavored-mass terms may be renormalized indepen- 
dently, and the (1,15) splitting could be broken. The calculation in |Q] shows that it is not the case, 
and the species-splitting is kept also in the interacting theory. 

2.2 Staggered- Wilson fermion 

There are also non-trivial flavored-mass terms for staggered fermions, which satisfy 75 her- 
miticity and work as a covariant laplacian. The staggered fermions with flavored-mass terms was 
referred to in the classical paper by Golterman and Smit and recently revisited by Adams in 
[|I|] as extension of Wilson and overlap fermions. This setup has been intensively studied in several 
works [||, ^, [7|, |9], [HJ. One advantage of this formulation is reduction of matrix sizes in the quark 
propagator, which would lead to reduction of numerical costs in lattice QCD with overlap fermions. 
In this subsection we discuss details of progress and problems in this course. 

For staggered fermions we only have two possibilities of flavored-mass terms satisfying e ~ 
75 <8> ^5-hermiticity: 1 <8> £5 and 1 ® a^ v . The flavored-mass term proposed by Adams has a former 
form, 

= e £ T7 1 T7 2 T73T7 4 C 1 C 2 C 3 C4 = (1 <8> &) + 0(a), (2.8) 

sym 

with {ri^xy = (-l)* 1+ " +r "- 1 <5 VjV and (e)^ = (-l) Xl+ ~ +M S x>y . where the factor 1/24 is hidden in 
the symmetric sum Y,sym. ■ The four tastes fall into £5 = + 1 subspace with two tastes and ^5 = — 1 
subspace with the other two tastes. Indeed, as shown in ((Tj, §], the degenerate Dirac spectrum is 
split into two branches. By introducing proper mass shift, the staggered version of Wilson fermion, 
which we call "staggered- Wilson fermion", is written as 

Srf = Y.XxiD^xyXy = E&fa/i^/i + K 1 J rM s s)+mo] xy Xy, (2-9) 

xy xy 



We re-derive (2.8) from the flavored-mass terms for naive fermions. As shown in \S\, spin diago- 



nalization decomposes Mp into four equivalent staggered flavored masses as 

^AQAQV* -> ±Zv(£T7imi3l4CiC2C 3 C4)Xv (2.10) 

Indeed, the Dirac spectrum of the staggered- Wilson operator shown in [|l|, ||] is the same as the left 
figure of Fig. ^| up to the degeneracy of the eige values. 

We now investigate symmetry of the staggered- Wilson fermion. We begin with review of the 
original staggered symmetry which is given by 

{Co.S^,/,,/?^} x {U e (l)} (2.11) 
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where Co is lattice charge conjugation, is is shift transformation, I s is spatial inversion, R^ v 
is hypercubic rotation, and U e {\) is residual chiral symmetry: % x — > e 6e *Xx- As well-known, the 
flavor and rotation symmetries are summarized as T4 xi SW^diag- I"4 is a Clifford group operating 
as flavor reflection and SW4.diag is the diagonal hypercubic subgroup of euclidian rotation SO(4) 
and flavor SU(4). By looking into its timeslice group it is shown that 15 pseudoscalars of 
fiavor-SU (4) fall into 7 irreducible representations (irreps) as 1 : §4,^45,^5, 3 : ^,^5,^,^4. Now 
let us move to the symmetry of the staggered- Wilson, which was intensively studied in [||, Bj. The 



flavored-mass (2.8) breaks the staggered symmetry down to 

{Co,^,R^v}, (2.12) 

where E'^ = E^T^. We take a brief look at this symmetry in terms of physical parity, charge con- 
jugation and the hypercubic symmetry. (See [jH, for details.) Regarding the parity, the action 
is invariant under the transformation ~ (74 <8> 1), which essentially stands for parity transfor- 
mation. For charge conjugation, the staggered charge conjugation symmetry Co remains intact in 
this case. Thus, physical charge conjugation for the two-flavor branch can be formed in a similar 
way to usual staggered fermions as shown in [§]. For Lorentz symmetry, the action is invariant 
under the staggered rotation R^ v and the shifted-axis reversal E' . These two groups are enough 
to form the proper hypercubic group SW4 as with the staggered fermion. The staggered- Wilson 
action thus possesses charge conjugation, parity and hypercubic symmetry. This result means that 
the staggered-Wilson satisfies the minimum conditions for practical use in lattice QCD. 

We now need to figure out advantages of staggered-Wilson and staggered-overlap over the 
originals. The staggered-Wilson fermion has been studied from several viewpoints, including the 
index theorem, Aoki phase, numerical costs and taste-breaking. 

Index theorem and overlap kernel : The study on the index theorem in the staggered-Wilson 
fermion via the spectral flow of the hernitian Dirac operator, which means the net number of near- 
origin eigenvalue crossings counted with signs of slopes, was first done by Adams in [[T]], and 
followed by [§, g |9|]. In this case the index of the Dirac operator is defined as Index (D -c /) = 
—Spectral flow(^/) with H^iin^) = eD^(mo), and it was shown that the index correctly detects 
the gauge topology in this setup. The stability of separation of low and high crossings in eigenvalue 
flow was also discussed, which indicates applicability of //^(mo < 0) as overlap kernel. In [^] the 
index theorem for the generalized Wilson including the Brillouin fermion was also studied. 

Aoki phase : The parity phase structure for the staggered-Wilson fermion was studied by 
using the Gross-Neveu model in [^] and strong-coupling lattice QCD [Bj. The GN study shows 
that the parity-broken phase and the 2nd-order critical line exist as with the case of the usual 
Wilson fermion, which implies possible Aoki scenario in this formulation. The strong-coupling 
lattice QCD also indicates the existence of the massless pion and the PCAC relation near the parity 
phase boundary. We expect the lattice chiral perturbation works to study this topic further. 

Numerical costs : One of possible advantages of the staggered-Wilson fermion is that it could 
reduce the computation expense of overlap fermions. The computational costs in the staggered- 
overlap quark propagator were well studied in According to this work, there is competition 
between an advantage from the 4-times smaller matrix to invert and a disadvantage from 4-link 
hopping terms in the flavored mass. Since the smaller matrix size requires fewer matrix-vector mul- 
tiplications in calculation of the sign function, the total CPU cost becomes smaller in the staggered- 
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overlap at least for the free theory. However, this merit is ruined by large gauge fluctuation due to 
the 4-link transporters. The gauge fluctuation reduces the splitting of the two branches, and unitary 
projection and inversion get more difficult. As a result, the reduction of the total CPU time is about 
0(2) for [5=6. Thus, so far, the numerical advantage of the staggered-overlap over the usual over- 
lap is not so significant. We however note that we can improve staggered-Wilson setups in terms 



of both staggered [ |19| ] and Wilson. The numerical cost would be reduced through the improved 
versions of the staggered-Wilson. 

Taste breaking : It is obvious that the two flavors in staggered-Wilson suffer SU (2) taste 
symmetry breaking. But it does not necessarily mean the mass-splitting of the associated three 
pions 7to,7i±. The pion mass splitting depends on the remaining discrete flavor symmetry: If 
this symmetry is large enough, we have degenerate three pions even at the finite lattice spacing. 
Recently Ref. [Jsj] has reported that the classification of operators by the timeslice symmetry indi- 
cates the three-degenerate pions: For pseudoscalar mesonic operators in staggered-Wilson, the 7 
irreps and flavor-singlet operators of the original staggered are mixed in <§ 5 pairs. And by focus- 
ing on the physical two-flavor branch £ we find the three pion states 7lo, 7i± are in the 3d irreps 
as i(j5 ® 0\)t (i = 1,2,3). The same work also shows that, in the pion potential of the stggered- 
Wilson chiral Lagrangian, the taste breaking starts from 0(a ) and 0(a 2 m). This result implis that 
the two flavors in the staggered-Wilson fermion suffer relatively small flavor breaking and suits the 
two-flavor lattice QCD. 

Another type: In [Q, the other type of staggered-Wilson fermions were proposed. This 
type corresponds to the case with 1 (g> a^ v , the only other possibility satisfying the staggered 75- 
hermiticity. This flavored mass is given by 

= /(T] 12 Ci2 + T734C34) = [l®(<Tl2 + 0-34)] + ^(fl), (2.13) 

with (T] M vW = ^vVtiVv8 x , y , (e MV )^ = (-1)^+^5^, C pv = (C M C V + C v C ll )/2. We note that 
the original flavored mass proposed by Hoelbling in [Q] has a different form, but we here discuss 
the above form since it has larger discrete symmetry [^, |[|. This flavored mass splits four tastes 
into three branches: one with positive mass, two with zero mass and the other one with negative 
mass. It thus produces single-flavor staggered-Wilson fermion, which seems to have wider appli- 
cability. However, the symmetry of the staggered fermion with this flavored mass term is small 
as {C7-,S^,/?i2,/?34,/?24^3i} with Cj = R2iR 2 l3 Co. Breaking of rotation symmetry indicates that 
we need to tune parameters to restore Lorentz symmetry. Indeed the recent study on symmetries 
of staggered-Wilson fermions by Sharpe [§| reports that recovery of Lorentz symmetry requires 
fine-tuning of parameters in the gluonic sector in lattice QCD. 



3. Central Branch 

In this section we discuss another way of use of flavored-mass or Wilson fermions. We usually 
concentrate on the one-flavor edge branch of the Wilson Dirac spectrum. However, the flavored- 
mass fermions without onsite terms (or species-singlet mass term) have a larger symmetry in gen- 



eral. In terms of the flavor-chiral symmetry of naive fermions in (2.2), the onsite term (~ y n y n ) 



breaks the invariance under any transformation of the axial-type group F x in Eq.(2.4). Thus, 
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dropping onsite terms can restore some invariance under the group, and the action comes to have a 
larger symmetry, including axial symmetry. 

In [0] and JToj] , it was shown that the Wilson fermion with the condition rriQ + Ar = has an 
extra U(l) symmetry besides the usual U(l) baryon symmetry. The condition corresponds to the 
six-flavor central branch of the Wilson Dirac spectrum as shown in Fig. |3| We call it "central-branch 
fermion". The action in this case is given by 

Scb = £ YnY^iD^Yn ~ r£ YnCuVn- (3.1) 

n,H n.jl 



In this case the U(4) x U (4) symmetry of naive fermions in ( |2.3| )(|2.4[) are broken down to 

r« -»■ {14, rH -> {(_i)«i+-+»4}. ( 3.2) 

The usual Wilson fermion has only the vector symmetry under I4 £ T^ + \ The invariance under 
(— 1)" 1+ - + " 4 € r(> is restored only with the central-branch condition niQ + Ar = 0. It is notable 
that the form (— i)»i+-+»4 j s the same as the staggered chiral transformation U e {\), which works 



as the axial rotation as shown in ( |2.2| ). It is obvious that the quark mass term is not invariant 
under this, and the extra symmetry prevents the quark mass term from being generated via 
loop effects. It means that the additive mass renormalization cannot occur in this case. Although 
the central branch of the usual Wilson fermion has six flavors, the symmetry enhancement on the 
central branch is generic with the flavored-mass fermions and this fact inspires us to consider new 
realizations of chiral fermions. 

The one-loop lattice perturbation in [^] shows that the additive quark mass renormalization 
becomes zero unlike the usual Wilson. The quark self-energy at one-loop level in a massless case 
is in general written as 

80 /% 
167T 2 V a 

(cc) 

We now focus on Zo, which is composed of the sunset L (sun) and tadpole (tad) diagrams. 
a = 1,2,..., 6 identifies the six poles of the propagator, which we denote as %y = (0,0,7T,7r), 
ftp = (0, 71,0, 7l), ■ ■ ■. For the central-branch fermion (fO|), Z^^sun) is calculated as 



+ iJnP^i • (3-3) 



(«)/ ^_sl r [ n d 4 k ^ (sin 2 ^-cos 2 ^)(-£ A cos£;0 + sin 2 fc p 



zW(sun) = ^C F / — y 1— 2 — ^Li^ 1 - 2E e <> = , (3.4) 
^-( 27r ) 4 r4(LAsin 2 ^)(E,sm 2 ^ + (-L,cos^) 2 ) 

• («) 

where e m p takes ±1 depending on or % locations of the poles in the direction p. This sign 
difference leads to cancelation between dimensions p in the integral for any pole a. The total 
contribution from the tadpole diagram, which we denote /( a )(tad), is given by 

' (a) M)= f fX l^lTt £ M*Pp + 4 K) ) + ^s(ap p + 4*>j) oc i7ilPll 
J-k{27i) 4 8a£ A sin / ^ A /2 y V p ) 

(3.5) 

(cc) ( cc) 

where we again have cancelation between dimensions p in £ () (tad), which means L (tad) = 

for any of the six poles. In the end the total £0 becomes zero for each of the six Dirac poles as 

= Eq (sun) + L^\ta.d) = 0. This verifies absense of additive renormalization in the setup. 
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Figure 4: The diagrams contribut- 



Figure 3: The Wilson Dirac spec . 

. „ r, t,, ( ,ing to the 1-loop fermion self- Figure 5: The Aoki phase diagram 
trum with mo + 4r = 0. The central 

branch crosses the origin. 



energy. The additive mass renor- with Wilson fermions. The red line 
malization from these diagrams be- corresponds to the central branch, 
comes zero. where the extra symmetry emerges. 



One concern on this setup is how the parity-broken phase affects the central branch. The 
central branch (mo + 4r = 0) corresponds to the central cusp in the parity phase diagram in Fig. [|. 
It is believed that this point is located within the parity broken phase at strong and middle gauge 
couplings, and the six flavors may be mutilated. Not to mention, we do not need to consider this 
problem if we take the chiral limit from the parity-symmetric phase to the central branch as shown 
in Nevertheless, we will find an important aspect of the central branch fermion by looking into 
the relation with the Aoki phase. The strong-coupling and large-A^ limit of lattice QCD right on 



the Wilson central branch was studied in QIOJ. It shows that the condensates in the strong-coupling 
limit of the central branch are given by 

_ M w _ 1 

°~4r 1, 71 ~ 16^(1 +r 2 



.(8r 4 -A4(l + ^)), 



(3.6) 



where a and % essentially stand for chiral ( and pion ( YYs Y) condensates with Mw = mo + Ar. 
It also shows that one of mesonic excitations has mass in the form as 



cosh (mspA 



1 + 



2Ml{\6+Ml) 



16 



15M 2 , 



(3.7) 



For the central branch Mw = m^+Ar = 0, we have a = 0, n = 1/2(1 +r 2 ) and mspA = 0. This result 
means that the extra symmetry is spontaneously broken and it leads to a massless NG meson, which 
seems to mimic QCD successfully. However, the symmetry is not broken by the chiral condensate 
(Yy) but by the parity broken condensate (YJsY)- Rather, the chiral condensate becomes exactly 
zero on the central branch. To summarize, Strong-coupling lattice QCD on the central-branch 
fermion has the condensate as (Yy) = 0, (YYsY) 7^ 0- This is a situation as if the roles of operators 
YY an d ^75^ are exchanged. In other words, the mass basis in this case is different from the usual 
QCD and we need to introduce a quark mass term in a twisted form, e.g. mYj5 Y- This fact indicates 
that the central-branch fermion is deeply related to the twisted-mass Wilson fermion. We note that 
the maximally-twisted-mass Wilson fermion is regarded as the average of the two one-flavor edge 
branches, and the central branch is located at the center point between these two. It implies that 
the central branch fermion can be seen as an automatic realization of the maximally-twisted-mass 
Wilson fermion without 0(a) discretization errors. If this view is correct, we can expect the parity 
breaking disappears in the continuum limit as in the twisted-mass Wilson. 
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We have another concern on this case : positivity of determinant. In the case of the maximally- 
twisted-mass Wilson QCD, the traceless flavor matrix T3 in the twisted-mass term ixj/Ys T3 y works to 
keep the quark determinant positive and avoid the theta term. Unlike this case, the quark condensate 
(YYsY) m tne strong-coupling limit of the central branch has no traceless structure in the 6-species 
flavor space. This problem is deeply related to the question how the U(l) problem is incorporated 
in this formulation. However, we note the chiral symmetry in this setup is related to the generator 
(— 1)" 1+ - + " 4 instead of naive 75. Thus, if the true condensate at the weak coupling has a form 
\j/(—l)' n+ - + " 4 Y an d it has traceless flavor structure for the six flavors, we may not necessarily 
have the issue. We have no consensus on this problem yet, and we hope detailed numerical study 
can answer this question. In any case, we can bypass the issue by introducing the pair of the central- 
branch fermions with T3 introduced in the mass term. In this case, the setup contains 12 flavors. 
This can be an alternative way for studying 12-flavor SU (3) gauge theory, whose existence of the 
infrared fixed point is now controversial. Not to mention, we are free from the problem when we 
take the chiral limit from the parity-symmetric phase as we discussed. 

So far we have concentrated on the central branch in the usual Wilson fermion. However, 
the flavored-mass fermions discussed in the previous section in general enjoy the extra symmetry 
restoration when the onsite term is dropped. For example, the fermion action with Mp but without 
mass shift has larger symmetry than (j2T6|): the flavor-chiral symmetry of \jf n (Y^D^ + Mp)y n is 

r (+) ^|i 4 , (-i)' l|+ - + " 4 7 5 , (-i) Vvi %^}, r(-)^{(-i)% T „ (-i)^y M y 5 }. (3.8) 

Restoration of the invariance under the axial-type group r'~' indicates absence of additive mass 
renormalization again. However, in this case, we have no central branch as seen in the left figure 
of Fig. || and we cannot take advantage of the larger symmetry. The staggered- Wilson fermion in 



(2.9) also has the similar property: in (2.9) without mass shift has larger symmetry as 



{Co, C r E M , C? T I„ R^}. (3.9) 

where C' T is given as a special charge conjugation C' T : % x — > %l \ Xx — > Xl > ~* ^uj Again, 
this case is not useful since there is no central branch. However the situation is different for the 
other type of the staggered-Wilson fermion in (2.13) [^, ^J: the symmetry of the action (2.13) 



without mass shift becomes larger as 

{Cr, C' T , S|j, /?i2, R34, ^24^31}- (3.10) 

The point is that we have a two-flavor central branch in this case. The extra symmetry is special 
charge conjugation C' T , which prohibits the additive mass renormalization. The numerical calcula- 
tion for this case was performed in [^]. Although this fermion has possibility of being a two-flavor 
setup without necessity of mass parameter fine-tuning, we have necessity of other parameter tun- 
ings for Lorentz symmetry for this case as we discussed [gj, which compensates the advantage. For 
our reference, the discrete symmetry of the staggered kinetic term, the two flavored-mass terms and 
the usual mass term are summarized in Table]]]. 

The interesting goal in this direction is to construct a two-flavor central-branch fermion by 



using a proper combination of the four flavored-mass terms shown in (2.5). In principle, it is 
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Table 1: 

possible to find a combination of flavored mass which has a two-flavor central branch. However, 
the question is whether hypercubic symmetry holds in such a case. If it is broken as with the single- 



flavor staggered- Wilson fermion in ( |2.13| ), it is not useful. We need further study to answer this 
question. 

4. Minimal-doubling 

In this section we discuss another formulation with the desirable number of flavors and chiral 
symmetry. One drawback in the flavored-mass fermions in Sec. |2| and || is that it explicitly breaks 
all the axial symmetry except on the central branches. The interesting question is whether we can 
reduce species to two with keeping the remnant of chiral symmetry. We now consider free Dirac 
operators with the flavored mass Mp multiplied by ij4, 

aD fc (p) = ijp sm Pfl a + iy 4 M F (/?), (4.1) 

where we assume Mp has a hypercubic or cubic symmetric form. Here the degeneracy of 16 species 
is not lifted by the flavored-mass term, but by specie-dependent imaginary chemical potential terms 
which we call "flavored-chemical-potential(FCP)" Jl6|]. Although we can consider a real type of 



FCP terms, we concentrate on the imaginary type to avoid the sign problem. (See [16, 17] for real 



FCP.) It is notable that this Dirac operator has ultra-locality and exact chiral symmetry, 

{y 5 ,O fc (/7)} = 0. (4.2) 

In this setup it is possible to reduce 16 species to two without losing all chiral symmetries, which 
is the minimal number allowed by the no-go theorem. On the other hand, the chemical potential 
term breaks the lattice discrete symmetry to the subgroup as cubic symmetry, CT and P [O], 



which give rise to necessity of fine-tuning several parameters for a correct continuum limit [15]. In 
other words, instead of chiral symmetry breaking in Wilson fermion leading to the mass parameter 
tuning, FCP fermions have breaking of discrete symmetry leading to necessity of other parameter 
tuning. To look into details of this setup, we introduce one concrete example 

r 3 

Skw = £ Vn7ii D nYn + moWnVn + r L WYa{ 1 ~ c jWn + ^WnmWn + d A \j/ x y A D A Yn ,(4-3) 

« L 7=1 

where we introduce a Wilson-like parameter r. The dimension-3 and dimension-4 counterterms 
with parameters jj,^ and J 4 are also introduced for the later discussion on tuning. The associated 
free and massless Dirac operator in momentum space is 

4 3 

aD KW (p) = i X y ll siaapn + iY4(ll3 + 3r — r V cosapj + ^sina/^). (4.4) 

ii=i j=\ 
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Figure 6: Species-splittings in Wilson and Karsten-Wilczek fermions. Circled numbers stand for the number 



of massless flavors on each point. [ 16 1 



For simplicity we first consider /I3 = 0, d.4 = and r = 1. In this case we have only two zeros of 
the Dirac operator p = (0,0,0,0) and p = (0,0,0, 71 /a) while the other 14 species have 0(\/a) 
imaginary chemical potential (~ ^YJYAYn) an d are decoupled in the naive continuum limit. In 
Fig. H we compare species-splitting of KW fermions in the chemical potential direction to that of 
Wilson fermion in the mass direction. For general values of /X3, the number of physical flavors 
depends on ^3 as shown in the left of Fig. |7|: flavor in /X3 > 1, 2 flavors in — 1 < fa < 1, 6 
flavors in — 3 < JU3 < — 1, 6 flavors in — 5 < /I3 < — 3 2 flavors in — 7 < /i 3 < —5, flavor in 
/I3 < —7. The action (4.3) is known as the Karsten-Wilczek (KW) fermion [Oj, which is the first 



type of "minimally doubled fermions" [11, O, 13]. Although we concentrate on the KW fermion 



all through this section, there are other types including Borici-Creutz [12] and Twisted-ordering 



[]13j fermions. The properties shown in this section are qualitatively common with all types. 

The symmetry of minimally doubled fermions was intensively studied in Ref.[[l4|]. As we 



have seen in a generic form of the action (|4JJ), the KW fermion has one exact chiral symmetry, 
ultra-locality, spatial cubic symmetry, CT and P. Because of this lower discrete symmetry, we need 
to fine-tune parameters to take a correct continuum limit. The process of the tuning was well 



discussed in [15]. In the case of KW fermion, it is pointed out that we need to tune parameters 



for a dimension-3 {yijuy) and two dimension-4 (xj/y^xj/, F^F^) counterterms. Two of these 



counterterms are shown in the action (1.3), and we need one more counterterm for the plaquette 
action. The necessity of three-parameter tuning seems to compensate the merits of the minimal- 
doubling fermions (chiral, ultra-local and two-flavor), and they have not been used widely. 



However, as seen from the generic argument below (|4.1|), minimal-doubling fermions are seen 
as lattice fermions with species-dependent chemical potential terms. It is notable that the discrete 
symmetry of minimal-doubling fermions is the same as the usual hypercubic-symmetric lattice 
fermions with the chemical potential, e.g. naive fermions with the chemical potential Sn(ju) = 
UniVnYjD jYn + YnY4 ^n,n+4 Wn+4 — ^ ^nfi-^n-A) /2), where the discrete symmetry is broken 
down to spatial cubic symmetry, CT and P again. Therefore, from the viewpoint of universality 
class, it is likely that both setups belong to the same universality class. We expect the symmetry 
of lattice QCD with minimal-doubling fermions will be enhanced to that of the two-flavor in- 
medium QCD (spatial rotation, P, CT, SU (2)flavor-chiral) in a continuum limit. Thus it is natural 
to regard the minimal-doubling fermion as the finite-density system and consider the application to 
the finite-(r,/») QCD. When an effective way of bypassing the sign problem is found in the near 
future, the two-flavor lattice fermion suitable for the finite-density system with ultra-locality and 
chiral symmetry will play an important role. 

Nevertheless, the minimal-doubling fermions and naive fermions with chemical potential have 
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Figure 7: (Left) The number of flavors in free Karsten-Wilczek fermion as a function of the parameter /I3. 
See fll6| ] for details. (Right) Conjecture on /13-g 2 chiral phase diagram for Karsten-Wilczek fermion with 
r=landd 4 = 0[||]. 



difference in the way of introducing the chemical potential. As seen from the discussion in the 
beginning of this section, the minimal-doubling fermion decouples 14 doublers with 0(1/ a) imag- 
inary chemical potential by introducing the FCP term, which is roughly given by the Wilson term 
ij/Mp y incorporating z%. Thus, instead of the additive mass renormalization in Wilson fermion, the 
unphysical 0(1/ a) chemical potential renormalization occurs in the interacting minimal-doubling 
fermion even if we consider it as the in-medium system. Then we need to fine-tune the dimension-3 
parameter /X3 to realize physical 0(1) chemical potential as we needed to fine-tune the mass pa- 
rameter mo to realize physical quark mass in Wilson fermion. This necessity of /X3 tuning is also 
understood from the well-known fact that the naive introduction of chemical potential on the lat- 
tice gives rise to the unphysical divergence of energy density and requires a counterterm since the 
violation of the abelian gauge invariance takes place [|]]] . 

This large renormalization leads to instability of the number of flavors in the interacting theory 
since the /I3 dependence of flavor number in the left of Fig. [7] in a free theory changes due to 
the renormalization. In Ref. [16], the jX^-g 2 chiral phase structure for KW fermion is studied 
by using the Gross-Neveu model and strong-coupling lattice QCD (the right figure of Fig. [7]). 
As seen from this, the /»3 range with physical flavors gets narrower as the gauge-coupling gets 
larger, and the minimal-doubling range also moves and gets narrower with nonzero gauge coupling. 
It means that we at least need to adjust /I3 within the minimal-doubling range to keep the two- 
flavor property irrespective of the way of use. On the other hand, it is not clear whether we need 
the dimension-4 parameter tuning in the application to in-medium QCD. As shown in [|16|], for 
general values of /I3 and J4, the dispersion relation in the Dirac operator is modified as D(p) ~ 



ijiPi + z/4P4 y (1 + oIa) 2 — /»| + 0(ap 2 ). It indicates that we may need to tune to correct the 
dispersion relation for given value of /»3 even if we consider the finite-density system. 

So far it is not clear how many parameters should be tuned in the application of minimal- 
doubling fermions to the in-medium QCD. However, Ref. [17] shows that the QCD phase diagram 
derived from this setup is consistent with the phenomenological predictions. In Ref. the 
strong-coupling limit of lattice QCD with KW fermion in the presence of temperature and density 
was studied. We depict the results in Figs. [8] and |9[ The (T, Hb) phase diagram in Fig. [8] is consis- 
tent with strong-coupling lattice QCD with staggered fermions, while there are some quantitative 
differences. For example, the ratio of the transition baryon chemical potential at T = to the crit- 
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m = 0, T = 0.3 m = 0, T = 0.2 




Figure 8: Phase diagram for the chiral transition with 

. «nj n £ j n jj Figure 9: Chiral condensate c and the baryon density 

r= l,u.-i = —0.9 and m = for d& = 0. Green and red " 11 

,. Z oa A-i tt t - v f 1 Pb for (left) T = 0.3 and (right) T = 0.2 with t/ 4 = 0. 

lines show 2nd and 1st transition lines, respectively. v ° 

™ , ... j • u t * .< Top and bottom panels show the massless m = case 

The transition order is changed from 2nd to 1 st at the r r 

...... . , „ tr ;\ / n „„, noiA \ Hi"=ni with 1st and 2nd phase transitions and massive m = 

tncntical point {jj,g,T ln j = (0.804,0.234) [|17[. F . — , 



0.1 case the crossover and 1 st-order transitions [|17| 

ical temperature at jX B = is given by R° = }X C {T = 0)/T c (hb = 0) ~ 2.2 for KW fermion while 
R^ t ~~ 1 for staggered. In the real world, this ratio is larger, R° ~ 5.5, and larger R° with KW 
fermion in the strong coupling limit may suggest smaller finite coupling corrections. The location 
of the tricritical point in KW fermion is given by /?^ w ~ 3.4 and while R^ =~ 2.0 for staggered 
fermion. It is consistent with the recent Monte-Carlo simulations, which implies that the critical 
point does not exist in the low baryon chemical potential region, Hb/T < 3. These observations 
reveal the usefulness of KW fermion for research on QCD phase diagram. 

In the end of this section, we comment on the question whether the other 14 flavors are re- 
ally decoupled in the continuum limit of minimal-doubling fermions. If fermions with infinitely 
large chemical potential are not decoupled, there should be some unphysical influence on the two- 
flavor QCD. This problem is discussed in [p2|], and it is pointed out that the subtraction of the free 
energy might be necessary to derive physical thermodynamical quantities. We still need further 
investigation in the application of this type of lattice fermions. 



5. Summary 

In this review paper we discuss the recent progress of lattice fermions and their possibilities 
of improving lattice simulations. The formulations not only open a new course of lattice study, but 
also give us feedback for the existing simulations. 

The flavored-mass terms, which are generalizations of the Wilson term, produce new types of 
Wilson, domain-wall and overlap fermions. The proper sum of the flavored-mass terms leads to 
the improved Wilson fermions including the Brillouin fermion. By introducing flavored mass into 
staggered fermions we obtain staggered- Wilson fermions, which have the possibility of reducing 
overlap numerical costs and staggered taste breaking. Although some numerical tests for these se- 
tups have not shown remarkable advantages, we expect the further tests with the improved versions 
could reveal their significance. The central branch in the flavored-mass fermions including Wilson 
fermion have the extra axial-type symmetiy. This symmetry implies no necessity of fine-tuning 
of the mass parameter in lattice QCD with this type. The lattice perturbation indicates that we 
have no additive mass renormalization in the interacting theory. From the strong-coupling lattice 
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QCD it is shown that this extra symmetry is spontaneously broken with the associated NG boson 
emerging. Although we need further investigation to judge its practical applicability, the interest- 
ing goal in this course is a two-flavor central-branch fermion, which can be a two-flavor setup with 
keeping chiral and hypercubic symmetry. Although the minimal-doubling fermion is a chiral and 
ultra-local setup, we need to tune the several parameters to restore Lorentz and C, P, T invariances 
in the continuum limit. This fermion can be seen as a lattice fermion with flavored chemical po- 
tential terms. This interpretation gives us new and nonperturbative understanding on the tuning 
of the dimension-3 parameter based on the chiral phase structure. Its use in the finite-density and 
-temperature systems was investigated and the strong-coupling QCD shows that it produces the 
phase diagram consistent to the phenomenological prediction. It requires numerical study to figure 
out requisite tuning process for this case. 
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